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We present a systematic study of the production of the heavy quarkonium, i.e., |(cc)[n]) 

, |(ftc)[n]) (or |(c6)[n]}), and |(6fo)[n]} quarkonium [|(Qp')[ n ]) quarkonium for short], through 
Z° boson semiexclusive decays with new parameters [l| for the heavy quarkonium under the 
framework of the nonrelativistic QCD, where [n] stands for n 1 So, n 3 S i, n 3 Po, n 3 Pj (n = 1, • • ■ ,6; 

J = [0,1,2]). “Improved trace technology” is adopted to derive the simplified analytic expressions 
at the amplitude level, which shall be useful for dealing with these decay channels. If all higher 
I (QQ')[ n ]) quarkonium states decay to the ground state |(Q<5 , )[l 1 >S'o]) with 100% efficiency 
via electromagnetic or hadronic interactions, we obtain V(Z° —> |(cc)[l 1 So])) = 1476 KeV, 

T(Z° ->• |(6c)[l 1 So])) = 1485 KeV, T(Z° ->■ |(6fo)[l 1 5 0 ]» = 127.5 KeV. At the LHC and ILC with 
the luminosity C, oc 10 34 cm _2 s^ 1 , sizable heavy quarkonium events can be produced through Z° 
boson decays; i.e., about 5.9 x 10 5 (cc), 6.0 x 10 5 (be) [or ( cb )], 5.1 x 10 4 (bb) events per year can 
be obtained. 

PACS numbers: 12.38.Bx, 13.66.Be, 14.40.Nd, 14.40.Pq 


I. INTRODUCTION 

The impressive volume of experimental data on multi¬ 
ple heavy quarkonium production were obtained by the 
experiments at the LHC in the last years. The LHCb and 
CMS Collaboration experiments have published studies 
of the B c -meson production and of the double < 7/\f r pro¬ 
duction 0 [S]. Also very interesting is that the candidate 
of the .B C (2S) state has been observed by ATLAS 0. 

At the LHC 0 , the International Linear Collider (ILC) 
00, and the newly purposed Z factory [s], the Z° bo¬ 
son is being and will be copiously produced. The Z° 
boson production cross section is about 34nb at the LHC 
with collider runs at the center-of-mass energy \/S = 14 
TeV and will be about 30nb at the ILC runs at the Z° 
pole energy 0 ], There will be about 10 9 Z° boson events 
being produced per year at the LHC and ILC with the 
luminosity C oc 10 34 c?rU 2 s -1 . Therefore, investigation 
of the heavy quarkonium production through Z° decays 
is worthwhile and meaningful. The study of the heavy 
quarkonium, e.g., |(cc)[n]}, |(&c)[n]) (or |(c6)[n])), and 
|(66)[n]) quarkonium, can help us to achieve a deeper un¬ 
derstanding of the quantum chromodynamics (QCD) in 
both the perturbative and nonperturbative sectors. The 
heavy quarkonium is presumed to be a nonrelativistic 
bound state of the heavy quark and antiquark. A won¬ 
derful theoretical tool to deal with the processes involv¬ 
ing quarkonium is the nonrelativistic quantum chromo¬ 
dynamics (NRQCD) 0], in which the low-energy interac¬ 
tions are organized by the expansion in v, where v stands 
for the typical relative velocity of the heavy quark and 
antiquark inside of the heavy quarkonium. The heavy 
quarkonium production itself is very useful for high pre- 
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cision physics in the electroweak sector and testing the 
perturbative QCD (pQCD) flo| - [l2l ]. For example, since 
its discovery by the CDF Collaboration [l3( , the B c me¬ 
son being the unique “doubly heavy-flavored” meson in 
the standard model (SM) has aroused great interest. The 
direct hadronic production of the B L meson has been 
studied systematically in Refs. [14i4l9| . For compensa¬ 
tion, it would be helpful to study its indirect production 
mechanisms. A systematic study of the B c meson pro¬ 
duction through the t quark or f quark, W±, and the 
Z° boson decay can be found in the literature IlL 120H27| . 
Because of sizable t/t quark, W ± and Z° boson events 
shall be produced at the LHC and ILC, the production of 
B c through their decay shall be helpful for determining 
the B c meson properties. 

Due to a high collision energy and high luminosity at 
the LHC and ILC, sizable amounts of the heavy quarko¬ 
nium events can be produced through Z° decays 0,0. 
So these channels may be an important supplement for 
other measurements at the LHC and ILC. Here we will 
further discuss the production of even higher \(QQ')[n]) 
quarkonium Fock states through the Z° semiexclu¬ 
sive decays with new parameters 0, where \(QQ')[n\) 
stands for |(QQ , )[n 1 S'o]), \{QQ')[n 3 S\]), |(QQ , )[n 1 Pi]), 
\(QQ')[n 3 Pj]) quarkonium [with n = 1, - - -, 6; J = 
(0,1,2)], respectively, and Q or Q' is c(c) or bfb) quark, 
respectively. As shown in Refs. [25, 0, HU, due to 
the color suppression of the amplitude and the relative- 
velocity suppression of the color-octet matrix element, 
the color-octet \(QQ')g[rC£) component provides negligi¬ 
ble contributions. We will only discuss the color-singlet 
state production channels. 

Under the NRQCD framework Q, a doubly heavy 
quarkonium is considered as an expansion of various Fock 
states. The relative importance among those infinite in¬ 
gredients is evaluated by the velocity scaling rule. To deal 
with heavy \(QQ'){n\) quarkonium production through 
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Z° semiexclusive decays, one needs to derive the pQCD 
calculable squared amplitudes. But the analytical ex¬ 
pression for the usual squared amplitude |£| 2 becomes 
too complex and lengthy for more (massive) particles in 
the final states and for higher-level Fock states to be gen¬ 
erated for the emergence of massive-fermion lines in the 
Feynman diagrams, especially to derive the amplitudes 
of the P- wave states. On account of needing to get the 
derivative of the amplitudes over the relative momentum 
of the constituent quarks for the P- wave states. It has 
been found that to do the numerical calculation using 
this conventional squared-amplitude technology becomes 
time-consuming for these complex processes, since the 
cross terms of the matrix elements increase with the in¬ 
crement of Feynman diagrams |£| 2 = S, :j MiM*. where 
i and j stand for the number of Feynman diagrams of 
the process. To solve this problem, the “improved trace 
technology” is suggested and developed in the literature 
[lL 123 - l27| ; it deals with the process directly at the ampli¬ 
tude level. After generating proper phase-space points, 
one first calculates the numerical value for the ampli¬ 
tudes, sums these values algebraically, and then squares 
the sum to get the squared amplitude, |£| 2 = |£iMi| 2 ; 
through this method, numerical simulation efficiency can 
be greatly improved in comparison to the conventional 
squared-amplitude technology. Moreover, under this ap¬ 
proach, many simplifications can be done at the ampli¬ 
tude level due to the fermion-line symmetries and the 
specific properties of each heavy quarkonium Fock state; 
then, one can even write down the analytic expressions 
for the amplitude. We will continue to adopt improved 
trace technology to derive the analytical expression for 
all the above-mentioned decay channels. 

The rest of the paper is organized as follows. In Sec. II, 
we introduce the calculation techniques for the Z° boson 
semiexclusive decays to \(QQ')[n\) quarkonium, where [n] 
stands for t^Sq, n 3 5i, n x Po, n 3 Pj (n = 1, - - -, 6; J = 
[0,1,2]). Then we calculate the production of |(cc)[?r]), 
|(6c)[n]), and |(66)[n]) quarkonium through Z° decay 
channels, i.e. Z° -A |(cc)[n]) + cc, Z° -A |(6c)[n]) + cb , 
and Z° -A |(&5)[n]) + bb , with new parameters [l| for 
the \(QQ')[n\) quarkonium in Sec. III. In Sec. IV, we 
discuss the production of |(QQ')[n]) quarkonium via Z° 
boson decays with the different number of flavor quarks 
rif under the five potential models. The final section is 
reserved for a summary. 


II. CALCULATION TECHNIQUES AND 
FORMULATION 

The four semiexclusive processes for the heavy quarko¬ 
nium production through Z° boson decays can be dealt 
with, i.e., Z° -A |(cc)[n]) + cc, Z° -A | (bc)[n\) + be (or 
Z° -A | (cb)[n)) + cb), and Z° -A |(66)[n]) +bb. These 
processes can be shorten for Z°(k) -a |(QQ')[ n ])(93) + 
Q(Q 2 ) + Q' (<?i), where Q or Q' stands for charm or beauty 
quark, and k and qi (* = 1,2,3) are the momenta of the 




FIG. 1: (color online). Feynman diagrams for the process 
Z°(k) -A \(QQ')[n])(q 3 ) + Q(q 2 ) + Q'(qi), where \{QQ')[n\) 
stands for |(QQ')K'S'o]>, KQQ'H^Si]}, ^QQ'^n 1 Pi]), 
\{QQ')[ r ^ Pj)) quarkonium (with n = 1, ■■■,6; J = (0,1,2)), 
respectively. Q or Q' is c(c) or b(b) quark, respectively. 

corresponding particles. According to the NRQCD fac¬ 
torization formula Q, the square of the semiexclusive 
amplitude can be written as the production of the per- 
turbatively calculable short-distance coefficients and the 
nonperturbative long-distance factors, the so-called non- 
perturbative NRQCD matrix elements. Its total decay 
widths dr can be factorized as 

dT = J2dT(Z° ^\(QQ')[n})+QQ')(0 H (n)), (1) 

n 

where (O h (n)) is the nonperturbative matrix element 
which describes the hadronization of a QQ' pair into the 
observable quark state H and is proportional to the tran¬ 
sition probability of the perturbative state QQ' into the 
bound state \{QQ')[n}). As for the color-singlet com¬ 
ponents, the nonperturbative matrix elements can be di¬ 
rectly related either to the wave functions at the origin or 
the first derivative of the wave functions at the origin [9], 
which can be computed via the potential models [1, 1291- 
IHj and/or potential NRQCD (pNRQCD) [lj, [34 l35| 
and/or lattice QCD [36) . 

The short-distance decay width T(Z°) can be ex¬ 
pressed as 

df(Z° -A l(QQ')N) + QQ') = (2) 

where means that one needs to average over the spin 
states of the initial particle and to sum over the color and 
spin of all the final particles. In the Z° rest frame, the 
three-particle phase space can be written as 

= <=>> 

The process to simplify the 1 —> 3 phase space with a 
massive quark/antiqark in the final state has been dealt 
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with in greater detail in Refs. [22], [24|. To shorten the pa¬ 
per, we shall not present it here, but the interested reader 
may turn to these references for the detailed technology. 
With the help of the formulas listed in Refs. [22j,[24|, one 
can not only derive the whole decay widths but also ob¬ 
tain the corresponding differential decay widths that are 
helpful for experimental studies, such as dT / dsi , dT/ ds 2 , 
dr/dcos 0 i 2 , and dr/dcos 0 i 3 , where Si = (<71 + q 2 ) 2 , 
S 2 = (< 7 i + ® 3 ) 2 , 612 is the angle between q\ and 92 , and 
013 is the angle between qi and < 73 . 

The hard-scattering amplitude for the specified pro¬ 
cesses can be dealt with as follows: 

Z° -A |(cc)[n]) + cc, Z° —> |( 6 c)[n]) + be , 

Z° -a \(bb)[n]) +bb. (4) 

The Feynman diagrams of the above three processes 
as Z°(k) -A | (QQ')( 93 ) + < 3 ( 92 ) + Q'(qi) is presented 
in Fig. [T] where the intermediate gluon should be hard 
enough to produce a bb pair or cc pair, so the amplitude 
is pQCD calculable. 

These amplitudes can be generally expressed as 

m 

iM =Cu si (q 2 )'^2A n v s 'j(qi), ( 5 ) 

n—1 

where m stands for the number of Feynman diagrams, 
s and s' are spin states, and i and j are color indices 
for the outing Q, Q' quark. The overall factor C stands 
for the specified quarkonium in the color-singlet, where 

C = 4 cosBw x °w stands for the Weinberg angle. 

A n in the formulas is similar in Ref. [25| , so do not listed 
here to shorten the paper. 

By using the improved trace technology, one can se¬ 
quentially obtain the squared amplitudes, and the nu¬ 
merical efficiency can also be greatly improved 0 , 112 - 
We adopt the improved trace technology to sim¬ 
plify the amplitudes M ss , at the amplitude level for the 
above-mentioned processes, the standard procedures for 
Z°(k) -A |(QQ')[n 1 )(( 73 ) + <3(92) + Q'{qi) are similarly 
presented in Ref. [251 ] . 

Finally, the decay widths over Si and S 2 can be ex¬ 
pressed as 

^ («> 

where mz is the mass of the Z° boson. The color- 
singlet nonperturbative matrix element (O h (n)) can be 
related either to the Schrodinger wave function ^(qq') (0) 
at the origin for the S -wave quarkonium states or the first 
derivative of the wave function V’(qq/)( 0 ) at the origin for 
the P -wave quarkonium states: 

(0 H (nS)) ~ |^|( Q q 0 [ „ s]> (0)| 2 , 


(o h (up)) ~ w im , )[nP]) m 2 - 

Since the spin-splitting effects are small, the difference 
between the wave function parameters for the spin-singlet 
and spin-triplet states at the same level are not distin¬ 
guished. The Schrodinger wave function at the origin 
'k|QQ , )[ n 5 ]} (0) and the first derivative of the Schrodinger 
wave function at the origin ^|(QQ,v n pi) (0) are related 
to the radial wave function at the origin R|(qq')[„ 5 ]>(0) 
and the first derivative of the radial wave function at the 
origin R'^qq^p]) (0), respectively, 

^ , |(QQ')[”S])(0) = \/l/47ri2|(QQ,)[ nS r])(0), 
^KQQ'Mn-P]^ 0 ) = a/ 3 / 47 ^ {QQ')[nP}) (°) ‘ ( 8 ) 

The radial wave function at the origin 7?|(QQ/)[ n 5])(0) 
and the first derivative of the radial wave function at the 
origin R^QQi^np-^ (0) relate to the number of active flavor 

quarks n/, the constituent quark mass of the |(Q<2 , )[ n ]) 
quarkonium, and the concrete potential models, respec¬ 
tively Q. Thus, i?|(QQq [nS ]}(0) and #[ (Q Q, )[nP]> (0) are 
adopted in Refs. [lj. 


III. NUMERICAL RESULTS 
A. Input parameters 

The input parameters are adopted as the following 
values flu Hg]: mz = 91.1876 GeV, Z° full width 
T{Z°) = 2.4952 GeV. 9w = arcsin -\/0.23119 is the Wein¬ 
berg angle, mw = 80.399GeV. We set the renormal¬ 
ization scale to be m( c5 ) and m^c) of |(cc)) and |(&c)) 
quarkonium for leading-order a s running , which leads to 
a s = 0.26 and of |(66)}-quarkonium for a s = 0.18. 

To ensure the gauge invariance of the hard amplitude, we 
set the \{QQ')[n)) quarkonium mass M to be toq + toq'. 
We adopt the values derived in Refs. 0,111 and list them 
in Table Q] since it is noted that the Buchmuller and 
Tye potential (B.T. potential) has the correct two-loop 
short-distance behavior in QCD and the decay 

widths are related to the constituent quark mass of the 
|(QQ')[n]> quarkonium. 


B. Heavy quarkonium production via Z° decays 

The decay widths for the \ (QQ')[n\) quarkonium states 
and the production channels through Z° decays, i.e., 
Z°_—> |(cc)[u]) + cc, Z° -A |(6c)[n]) + be ,and Z° —> 
| (bb)[n\) + bb, are listed in Tables HU IIII1 and IIVI within 
the B.T. potential 0 ]. 
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TABLE I: Mass of the constituent quark and radial wave functions at the origin under the B.T. potential pj. 


KQQ’M) 

Mass and wave functions 

n = 1 

n = 2 

n = 3 

n = 4 

n = 5 

n = 6 

n = 7 

S states of |(cc)[n]) 

m c {GeV) 

1.48 

1.82 

1.92 

2.02 

2.12 

2.25 



|P, [nS1> (0)| 2 (GeV 3 ) (n /= 3) 

2.458 

1.617 

0.969 

0.796 

0.701 

0.721 


P states of |(cc)[n]) 

m c {GeV) 

1.75 

1.96 

2.12 

2.26 

2.38 




l^[n P i)(0)| 2 (GeV 5 ) (n/=3) 

0.322 

0.224 

0.387 

0.467 

0.499 



S states of |(fec)[n]) 

m c {GeV) 

1.45 

1.82 

1.96 

2.10 

2.15 




m b {GeV) 

4.85 

5.03 

5.15 

5.30 

5.45 




|P| InS]> (0)| 2 (GeV 3 ) (n/=3) 

3.848 

1.987 

1.347 

1.279 

1.118 



P states of |(6c)[n]} 

m c {GeV) 

1.75 

1.96 

2.15 

2.26 





m b {GeV) 

4.93 

5.13 

5.25 

5.37 





l^[n P i)(0)| 2 (GeV 5 ) (n/=3) 

0.518 

0.500 

0.729 

0.823 




S states of |(fe&)[n]) 

m b {GeV) 

4.71 

5.01 

5.17 

5.27 

5.41 

5.50 

5.58 


|P| InS1> (0)| 2 (GeV 3 ) {n f = 4) 

16.12 

6.746 

2.172 

2.588 

2.665 

2.576 

2.377 

P states of |(bb)[n]} 

m b {GeV) 

4.94 

5.12 

5.20 

5.37 

5.47 

5.56 



\Ru nv]) m 2 {GeV 5 ) {n f = 4) 

5.874 

2.827 

2.578 

3.217 

3.573 

3.669 



TABLE II: Decay widths (in unit KeV) for the production of |(cc)[n]) (rif = 3) through Z 0 decays. 



n = 1 

n = 2 

n = 3 

n = 4 

n = 5 

n = 6 

Z° (cc)[n 1 5o] + cc 

313.0 

108.0 

54.67 

38.25 

28.90 

24.59 

Z° —> (cc)^ 3 ^!] + cc 

325.0 

112.3 

56.93 

39.86 

30.13 

25.66 

Z° —> (cc)[n 1 Pi] + cc 

27.84 

11.66 

11.91 

10.50 

8.473 


Z° —t (cc)[n 3 Po] + cc 

39.15 

14.81 

17.43 

15.12 

12.49 


Z° —> (ec)[n 3 Pi] + cc 

41.52 

16.22 

18.78 

16.34 

13.40 


Z° — t (cc)[n 3 P 2 ] + cc 

18.32 

5.857 

7.173 

6.303 

5.380 



TABLE III: Decay widths (in unit KeV) for the production of |(6c)[n]} (n/ = 3) through Z° decays. 



n = 1 

n = 2 

n = 3 

n = 4 

n = 5 

Z° (fee) [n 1 So] +bc 

342.2 

87.71 

47.22 

36.11 

29.22 

Z° (fee) [n 3 Si] + fee 

492.4 

119.5 

63.46 

48.02 

38.96 

Z° ->• (fec^Pi] + bc 

17.09 

9.237 

8.617 

7.592 


Z° ->• (fec)[n 3 P 0 ] + be 

11.69 

6.962 

6.954 

6.313 


Z° {bc){n 3 Pi]+bc 

21.41 

12.05 

11.41 

10.14 


Z° ->• (6c)[n 3 P 2 ] + bc 

20.97 

11.20 

9.954 

8.651 



TABLE IV: Decay widths (in unit KeV) for the production of | (66) [n]) (n/ =4) through Z° decays. 



n = 1 

n — 2 

n = 3 

n = 4 

n = 5 

n = 6 

n = 7 

Z° (6fe)[n 1 S 0 ] + fefe 

28.48 

9.537 

2.737 

3.039 

2.839 

2.581 

2.256 

Z° (feb)[n 3 Si] + bfe 

31.83 

10.74 

3.095 

3.446 

3.231 

2.944 

2.579 

Z° (fefe) [n 1 Pi ] + fefe 

1.387 

0.550 

0.461 

0.483 

0.484 

0.455 


Z° -)• (fefe) [n 3 P 0 ] + fefe 

2.199 

0.878 

0.739 

0.779 

0.785 

0.740 


Z° ->■ (fefe) [n 3 Pi ] + fefe 

2.104 

0.835 

0.700 

0.734 

0.738 

0.693 


Z° (fefe) [n 3 P 2 ] + fefe 

0.870 

0.345 

0.290 

0.304 

0.305 

0.287 
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FIG. 2: (color online). Differential decay widths dT/dsi (Left) and dV/ds 2 (Right) for Z° —► |(cc)[n]) + cc, where the diamond 
line, the dotted line, the solid line, the dash-dotted line, and the crossed line are for |(cc)[n]), where n = 1, ■ ■ ■, 5, respectively. 
|(cc)[n]) stands for the summed decay width of [nl 1 5o], [n 3 Si], [n 1 Pi], and [n 3 Pj] (J = 0,1,2) at the same nth level. 




FIG. 3: (color online). Differential decay widths dF/dsi (Left) and dV/ds 2 (Right) for Z° —¥ |(6c)[n]) + 6c, where the diamond 
line, the dotted line, the solid line, and the dash-dotted line are for |(6c)[n]), where n = 1, • • ■ ,4, respectively. |(6c)[n]) stands 
for the summed decay width of [n 1 So], [n 3 5i], [^Pi], and [n 3 Pj] (J = 0,1,2) at the same nth level. 


From Tables ITTlIIVl it is noted that, in addition to the 
ground 15-level states, the higher |(QQ')[ n ]) quarkonium 
states can also provide sizable contributions to the total 
decay widths. For convenience, we have used [n5] to 

I 


• For |(cc)[n]) quarkonium production through the 
decay channels Z° —> |(cc)[n]) + cc, the total de¬ 
cay widths for all 25, 35, 45, 55, 65, IP, 2P, 3P, 
4P, and 5P-wave states is 34.5%, 17.5%, 12.2%, 
9.3%, 7.9%, 19.9%, 7.6%, 8.7%, 5.9%, and 6.2% of 
the summed decay widths of r] c and J/'Sf. Con¬ 
sidering that the LHC runs at the center-of-mass 


present the summed decay widths of [n 1 5o] and [n 3 5i] 
at the same nth level, and [nP] to represent the summed 
decay widths of [n 1 Pi] and [n 3 Pj](J = 0,1,2) at the 
same nth level. 


energy y/S = 14 TeV and the ILC runs at the Z° 
pole energy with the luminosity C oc 10 34 cm _2 s _1 , 
one expects that about 1.0 x 10 9 Z° events per year 
can be generated. Then we can estimate the char- 
monium events generated through Z° decays; i.e., 
2.6 x 10 5 | (cc)[15]), 8.9 x 10 4 |(cc)[25]), 4.5 x 10 4 
|(cc)[35]), 3.2 x 10 4 |(cc)[45]), 2.4 x 10 4 |(cc)[55]), 
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FIG. 4: (color online). Differential decay widths dT/dsi (Left) and dT/ds 2 (Right) for Z° —> |(6&)[n]} +bb, where the diamond 
line, the dotted line, the solid line, the dash-dotted line, and the crossed line are for |(&6)[n]), where n — 1, ■ ■ ■, 5, respectively. 
|(&5)[n]) stands for the summed decay width of [n 1 i So], [n 3 5i], [n 1 Pi], and [n 3 Pj] (J = 0,1, 2) at the same nth level. 




FIG. 5: (color online). Differential decay widths dT / dcos0\2 (Left) and dT / dcosdis (Right) for Z° —> |(cc)[n]} + cc, where the 
diamond line, the dotted line, the solid line, the dash-dotted line, and the crossed line are for |(cc)[n]), where n = 1, • • •, 5, 
respectively. |(cc)[n]} stands for the summed decay width of [n 1 ^], [n 3 Si], [n 1 .Pi], and [n 3 Pj] (J = 0,1,2) at the same nth 
level. 




FIG. 6: (color online). Differential decay widths dY/dcos9\2 (Left) and dT/dcosOiz (Right) for Z° —> |(&c)[n]) +bc, where the 
diamond line, the dotted line, the solid line, and the dash-dotted line are for |( 6 c)[n]}, where n = 1, • • •, 4, respectively. | (foe) [n]) 
stands for the summed decay width of [n 1 S'o], [n 3 ^!], [n 1 Pi], and [n 3 Pj] (J = 0, 1,2) at the same nth level. 
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FIG. 7: (color online). Differential decay widths dT/dcosdi 2 
(Up) and dF/dcosOis (Down) for Z° — > |( 66 )[n]) +bb, where 
the diamond line, the dotted line, the solid line, the dash- 
dotted line, and the crossed line are for |( 66 )[n]), where n = 
1, • • •, 5, respectively. |( 66 )[n]) stands for the summed decay 
width of [n 1 So], [n 3 5i], [n 1 Pi], and [ n 3 Pj] (J = 0,1, 2) at the 
same nth level. 


2.0 x 10 4 |(cc)[65]), 5.6 x 10 4 |(cc)[lP]), 2.0 x 10 4 
|(cc)[2P]), 2.2 x 10 4 |(cc)[3P]), 1.5 x 10 4 |(cc)[4P]), 
and 1.6 x 10 4 | (cc) [5P]) quarkonium events per year 
can be obtained. 

• For |( 6 c)[n]) quarkonium production via Z° boson 
semiexclusive decays, the total decay widths for all 
25, 35, 45, 55, IP, 2P, 3P, and 4P-wave states 
is 24.8%, 13.3%, 10.1%, 8.2%, 8.5%, 4.7%, 4.4%, 
and 3.9% of the summed decay widths of B c and 
B* . At the LHC and ILC, 3.4 x 10 5 |(6c)[15]), 
8.4 x 10 4 | (6c)[25]), 4.5 x 10 4 |(6c)[35]), 3.4 x 10 4 
| ( 6 c)[45]), 2.6 x 10 4 | ( 6 c)[55]), 2.9 x 10 4 |( 6 c)[lP]), 
1.6 x 10 4 |(6c)[2P]), 1.5 x 10 4 |(6c)[3P]), and 1.3 x 
10 4 |(6c)[4P]) quarkonium events per year can be 
obtained. 

• For |( 66 )[n]) quarkonium production via Z° boson 
semiexclusive decays, the total decay widths for all 
25, 35, 45, 55, 65, 75, IP, 2P, 3P, 4P, 5P, 
and 6 P wave states are about 34.4%, 9.7%, 10.8% 
10.1%, 9.2%, 8.0%, 10.9%, 4.3%, 3.6%, 3.8%, 3.8%, 
and 3.6% of the summed decay widths of r)t and T. 


At the LHC and ILC, bout 2.4 x 10 4 |(66)[15]), 
8.1_x 10 3 | (66)[25]), 2.3 x 10 3 |(66)[35]), 2.5 x 10 3 
| ( 66 )[45]), 2.4_x 10 3 |(66)[55]), 2.2 x 10 3 |(66)[65]>, 
1.9 x 10 3 | ( 66 ) [75]), and summed up, 7.1 x 10 3 
| ( 66 ) [P]) quarkonium events per year can be ob¬ 
tained. 

To better illustrate the relative importance of different 
production channels, we present the differential distri¬ 
butions dT/dsi , dT/ds 2 , dT / dcos 0 \ 2 , and dT/dcosOis for 
the above-mentioned three processes in Figs. [2] [3] HI EJ El 
[7] We have used \ (QQ')[n]) to represent the summed de¬ 
cay width of (QQ')K5 0 ]), (QQ'^S!}), {QQ')[n}P^), 
and (QQ')[n 3 Pj]) (J = 0,1,2) at the same nth level in 
the figures, where n = 1, • • •, 5. These figures show ex¬ 
plicitly that the excited Fock states, i.e., excited states 
of \(QQ')[n]) quarkioum (n = 2 , • • •, 5), can provide siz¬ 
able contributions in comparison to the lower Fock state 
|(QQ , )[l 1 5o]) or |(QQ , )[1 3 5 i]) in almost the entire kine- 
matical region. 

As all of the higher excited heavy quarkonium nS and 
nP Fock states almost decay to the ground spin-singlet 5 
wave state KQQOtl^So]) via electromagnetic or hadronic 
interactions, we obtain the total decay width of Z° boson 
decay channels within the B.T. potential: 

r(Z° |(cc)[l 1 5 0 ]) +cc) = 1476 KeV, (9) 

T(Z° |(6c)[l 4 5 0 ]) + be) = 1485 KeV, (10) 

F(Z° -> |(66)[l 4 5o]) +66) = 127.5 KeV. (11) 

Running at the center-of-mass energy y/S = 14 TeV 

at the LHC [f| and at the Z° pole energy at the ILC 
6 . 0 ] with luminosity 10 34 cm _ 2 s _1 , one may expect to 
produce about 1.0 x 10 9 Z° boson per year. Then we 
can estimate the event number of \{QQ')) quarkonium 
production through Z° boson decays, i.e., about 5.9 x 
10 5 |(cc)[n]) quarkonium events, 6.0 x 10 5 |( 6 c)[n]) (or 
|(c 6 )[n])) quarkonium events, 5.1 x 10 4 |( 66 )[n]) quarko¬ 
nium events per year. Bearing in mind the situation 
pointed out here and the possible upgrade for the LHC 
(SLHC, DLHC, etc. [40[), and the newly purposed Z 
factory with luminosity 10 36 cto _ 2 s _1 f||, the possibil¬ 
ity to study \(QQ')[n}) quarkonium via Z° boson decays 
is worth serious consideration. 


C. Decay widths under five potential models 

In this subsection, we discuss the uncertainties for the 
\{QQ r )[n]) quarkonium production through Z° boson de¬ 
cays. For the present calculation, their main uncertainty 
sources include the nonperturbative bound state matrix 
elements, the renormalization scale hr, and the con¬ 
stituent quark masses raj, and m c . These parameters 
are the main uncertainty source for estimating heavy 
\{QQ’) [n])-qu£irk on ium production. In the present, we 
discuss the decay widths of \(QQ')[n}) quarkonium pro¬ 
duction through Z° boson decays under the five poten¬ 
tial models , i.e., the B.T. potential PI H3, the John 







TABLE V: Decay widths (in unit keV) for the |(cc)[n]} quarkonium production channel Z Q —>• |(cc)[n]) + cc, whose bound-state 
parameters are adopted in the five potential models in Ref. pj. 



B.T. nf=3 / 4 [30] 

J. nf=3 / 4 [41] 

I.O. nf=3 / 4 [42J 

C.K. nf=3 / 4 [33] 

Cor. [29] 

M = I'.Sb; 

313.0 / 298.5 

142.5 / 127.0 

71.95 / 76.28 

92.45 / 101.2 

124.2 

M = [l 3 5i] 

325.0 / 309.9 

148.0 / 131.8 

74.71 / 79.20 

95.99 / 105.1 

128.9 

W = pdsy 

108.0 / 90.83 

70.60 / 60.78 

36.67 / 38.07 

41.01 / 43.55 

59.38 

W = [2 3 5r] 

112.3 / 94.45 

73.41 / 63.20 

38.13 / 39.59 

42.64 / 45.28 

61.74 

W = f3 1 5 0 ] 

54.67 / 49.76 

55.57 / 47.17 

29.23 / 30.13 

31.48 / 32.95 

46.32 

W = [3 3 Si] 

56.93 / 51.82 

57.87 / 49.12 

30.43 / 31.37 

32.78 / 34.31 

48.23 

[n] = [4^0] 

38.25 / 38.11 

46.61 / 39.21 

24.65 / 25.32 

26.04 / 27.10 

38.78 

W = [4 3 ^r] 

39.86 / 39.71 

48.57 / 40.86 

25.69 / 26.39 

27.14 / 28.24 

40.41 

[n] = frSo] 

28.90 / 27.12 

37.76 / 31.61 

20.08 / 20.58 

20.93 / 21.67 

31.41 

W = [5 3 Si] 

30.13 / 28.27 

39.37 / 32.95 

20.93 / 21.46 

21.82 / 22.59 

32.75 

[n] = ^So] 

24.59 / 24.62 

34.38 / 28.68 

18.35 / 18.79 

19.00 / 19.58 

28.72 

M = [ 6 3 ^i] 

25.66 / 25.70 

35.87 / 29.93 

19.15 / 19.61 

19.82 / 20.43 

29.97 

N = [iR] 

126.8 / 129.6 

67.75 / 53.17 

20.88 / 22.45 

29.15 / 31.90 

35.84 

M = [2 P] 

48.55 / 49.85 

66.97 / 51.36 

21.46 / 22.54 

27.74 / 30.13 

36.63 

[n] = [3P] 

55.29 / 54.01 

62.44 / 47.43 

20.29 / 21.29 

25.29 / 27.29 

34.86 

[n] = [4P] 

48.26 / 48.99 

58.49 / 44.13 

19.22 / 20.15 

23.36 / 25.11 

33.07 

N = [5 P] 

39.74 / 40.94 

55.27 / 41.50 

18.40 / 19.11 

21.90 / 23.42 

29.95 

Sum. 

1476 / 1402 

1138 / 919.9 

510.2 / 532.3 

598.5 / 639.9 

841.2 


TABLE VI: Decay widths (in unit keV) for |(6c)[n]) quarkonium production channel Z° —> |(6c)[n]} + be, whose bound-state 
parameters are adopted in the five potential models in Ref. [l[. 



B.T. nf=3 / 4 [20] 

J. nf=3 / 4 [41] 

I.O. nf=3 / 4 [42] 

C.K. nf=3 / 4 [33] 

Cor. [29] 

M = [i 1 ^] 

342.2 / 356.5 

179.7 / 162.7 

552.3 / 467.9 

116.0 / 128.7 

158.6 

M = [i 3 Pi] 

492.4 / 513.0 

258.6 / 234.0 

794.8 / 673.3 

166.9 / 185.2 

228.2 

M = p 1 *] 

87.71 / 61.67 

79.68 / 69.17 

95.74 / 86.43 

46.17 / 49.22 

70.36 

M = [2 3 Pi] 

119.5 / 84.02 

108.6 / 94.24 

130.4 / 117.8 

62.91 / 67.06 

95.86 

[n] = [3 1 5o] 

47.22 / 42.38 

58.05 / 49.57 

45.61 / 41.58 

32.71 / 34.32 

51.32 

[n] = [3 3 Sr] 

63.46 / 56.96 

78.02 / 66.62 

61.29 / 55.88 

43.96 / 46.12 

68.97 

M = [4^0] 

36.11 / 36.56 

45.82 / 38.74 

26.57 / 24.42 

25.49 / 26.51 

40.71 

M = [4 3 <Si] 

48.02 / 48.62 

60.94 / 51.51 

35.33 / 32.48 

33.90 / 35.25 

54.14 

[n] = ^So] 

29.22 / 31.83 

41.06 / 34.47 

19.18 / 17.69 

22.69 / 23.44 

36.41 

[n] = [5 3 5r] 

38.96 / 42.44 

54.75 / 45.96 

25.58 / 23.59 

30.25 / 31.26 

48.54 

M = [IP] 

71.16 / 103.9 

56.74 / 45.47 

78.72 / 64.70 

25.55 / 28.71 

30.09 

M = [2P] 

39.45 / 34.40 

54.12 / 42.37 

38.11 / 32.35 

24.62 / 27.30 

29.98 

M = [3P1 

36.94 / 39.27 

47.78 / 36.94 

21.08 / 18.19 

19.76 / 21.58 

27.21 

M = [4P] 

32.70 / 36.91 

45.85 / 35.12 

14.46 / 12.59 

19.82 / 21.57 

26.54 

Sum. 

1485 / 1488 

1169 / 1007 

1939 / 1669 

671.5 / 726.2 

966.9 


L. Richardson potential (J. potential) [l, 41], the K. Igi 
and S. Ono potential (I.O. potential) fll. l42l. l43j . the 
Yu-Qi Chen and Yu-Ping Kuang potential(C.K. poten¬ 


tial) [l, Hj43) L and Coulomb-plus-linear potential (Cor. 
potential) [L, [29( in detail. 


The decay widths for \(QQ')[n]) quarkonium produc¬ 
tion under the five potential models with the different 


number of flavor quarks n / presented in Tables [V] IVI1 





















































9 


TABLE VII: Decay widths (in unit keV) for |(b6)[n])-quarkonium production channel Z° —> |(66)[n]) + bb, whose bound-state 
parameters are adopted in the five potential models in Ref. [l|. 



B.T. nf=4 / 5 [30] 

J. nf=4 / 5 [41] 

I.O. nf=4 / 5 [42J 

C.K. nf=4 / 5 [33J 

Cor. [29J 

M = i'.sV 

28.48 / 24.61 

12.51 / 9.827 

17.55 / 15.29 

9.314 / 10.69 

16.07 

M = [i 3 Si] 

31.83 / 27.51 

14.04 / 10.98 

19.61 / 17.09 

10.41 / 11.95 

17.96 

M = [2 1 S’o| 

9.537 / 10.49 

5.861 / 4.083 

4.894 / 4.262 

3.934 / 4.230 

6.745 

M = [2 3 Si] 

10.74 / 11.81 

6.601 / 4.598 

5.512 / 4.800 

4.431 / 4.763 

7.596 

M = [3 1 5o[ 

2.737 / 6.093 

4.286 / 2.845 

2.585 / 2.252 

2.797 / 2.930 

4.916 

M = [3 3 Sr] 

3.095 / 6.890 

4.846 / 3.218 

2.923 / 2.546 

3.163 / 3.313 

5.559 

M = [4 4 S 0 ] 

3.039 / 3.476 

3.578 / 2.314 

1.707 / 1.488 

2.316 / 2.392 

4.109 

M = [4 3 Si] 

3.446 / 3.941 

4.057 / 2.624 

1.936 / 1.687 

2.626 / 2.712 

4.659 

[n] = [5 4 5 0 ] 

2.839 / 2.377 

3.074 / 1.958 

1.218 / 1.063 

1.983 / 2.029 

3.541 

W = [5 3 5r] 

3.231 / 2.705 

3.499 / 2.228 

1.386 / 1.210 

2.256 / 2.310 

4.030 

[n] = ^So] 

2.581 / 2.251 

2.767 / 1.742 

0.943 / 0.824 

1.783 / 1.814 

3.189 

W = [6 3 <S'i] 

2.944 / 2.568 

3.157 / 1.987 

1.075 / 0.939 

2.034 / 2.069 

3.638 

N = [7 4 S 0 ] 

2.256 / 2.192 

2.540 / 1.585 

0.761 / 0.665 

1.636 / 1.656 

2.927 

W = [7 3 5r] 

2.579 / 2.506 

2.903 / 1.812 

0.870 / 0.761 

1.871 / 1.893 

3.346 

N = [IP] 

6.560 / 5.554 

1.836 / 0.986 

1.301 / 1.021 

1.241 / 1.501 

1.360 

M = [2 P] 

2.608 / 4.812 

1.980 / 0.987 

0.890 / 0.700 

1.221 / 1.427 

1.538 

[n] = [3P] 

2.190 / 3.411 

2.084 / 0.996 

0.675 / 0.531 

1.232 / 1.412 

1.666 

[n] = [4P] 

2.300 / 2.163 

2.031 / 0.946 

0.500 / 0.396 

1.170 / 1.326 

1.662 

[n] = [5P] 

2.312 / 2.053 

2.033 / 0.929 

0.402 / 0.319 

1.151 / 1.292 

1.691 

[n] = [6P] 

2.175 / 2.099 

2.034 / 0.915 

0.335 / 0.266 

1.135 / 1.266 

1.711 

Sum. 

127.5 / 148.4 

85.72 / 57.56 

67.07 / 58.11 

57.70 / 62.98 

97.91 


and IVIll where [nP] represent the summed decay widths 
of (QQ')[n l P\}) and ( QQ')[n 3 Pj ]} [with J = (0,1,2)] at 
the same nth level. The decay widths for the other four 
models are consistent with each other: taking the B.T. 
model decay width as the center value, for the channel 
Z° —> |(cc)[n]} + cc, we obtain the uncertainty for 

nf = 3 and (^@ 2 %) f° r n f = 4, where the lower value is 
for the I.O. potential; for the channel Z° —> \(bc)[n]) + bc, 
we obtain the uncertainty (-55%) for n f = 3 and (-51%) 
for rif = 4, where the upper value is for the I.O. poten¬ 
tial and the lower value is for the C.K. potential; for the 
channel Z° —> |(&&)[n]} + bb, we obtain the uncertainty 
(^ 55 %) f° r n f = 3 and (+ 5 ^%) for n/ = 4, where the lower 
value is for the C.K. potential. 

IV. CONCLUSIONS 

In this paper, we have made a detailed study on 
the |(cc)[n]}, |(&c)[n]), |( 66 )[?r])-quarkonium via Z° bo¬ 
son semiexclusive decays under the NRQCD framework, 
i.e., Z° —> |(cc)Jn]) + cc, Z° |( 6 c)[?i]) + be, and 
Z° —I \(bb)[n]) +bb, where [n] stands for [n 1 5o], [n 3 Si], 
[n 1 Pi], and [n 3 Pj], (n = 1, • • •, 7; J = [0,1, 2]). To pro¬ 
vide the analytical expressions as simply as possible, we 
have adopted the improved trace technology to derive 
Lorentz-invariant expressions for Z° boson decay pro¬ 
cesses at the amplitude level. Such a calculation technol¬ 


ogy will be very helpful for dealing with processes with 
massive spinors. 

Numerical results show that high excited states of 
\(QQ)[n}) quarkonium in addition to the ground IS wave 
states can also provide sizable contributions to heavy 
quarkonium production through Z° boson decays, so one 
needs to take the excited wave states into consideration 
for a sound estimation. If all the excited states decay to 
the ground state KQQ^^So]) with 100% efficiency via 
electromagnetic or hadronic interactions, we can obtain 
the total decay width for \(QQ')) quarkonium production 
through Z° boson decays as shown by Eqs. m-m - At 
the LHC at runs with center-of-mass energy \/~S = 14 
TeV and at the ILC at the Z° pole energy with the lumi¬ 
nosity £ oc 10 34 cm~ 2 s _1 , due to the high collision energy 
and high luminosity, sizable heavy quarkonium events 
can be produced through Z° boson decays; i.e., about 
5.9 x 10 5 (cc), 6.0 x 10 5 (be) [or(c&)[, 5.1 x 10 4 , (bb) 
events per year can be obtained. At the newly purposed 
Z factory with the high luminosity £ oc 10 36 cm _2 s _1 , 
the \(QQ')) quarkonium through Z° boson decays will 
be more abundantly produced. Therefore, one needs to 
take these higher excited states into consideration for a 
sound evaluation. 
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